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An ongoing limitation of terahertz spectroscopy is that the technique is generally limited to the study of 
relatively large samples of order 4 mm across due to the generally large size of the focal beam spot. We 
present a nested concentric parabolic reflector design which can reduce the terahertz focal spot size. This 
parabolic reflector design takes advantage of the feature that reflected rays experience a relative time delay 
which is the same for all paths. The increase in effective optical path for reflected light is equivalent to the 
aperture diameter itself. We have shown that the light throughput of an aperture of 2 mm can be increased 
by a factor 15 as compared to a regular aperture of the same size at low frequencies. This technique can 
potentially be used to reduce the focal spot size in terahertz spectroscopy and enable the study of smaller 
samples. 
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I. INTRODUCTION 

Time-domain terahertz spectroscopy (TDTSJP^ has 
seen tremendous growth is recent years and has found 
use in many different areas in clud ing security appli¬ 
cations and biohazard detectioiP^, detection of pro¬ 
tein conformational changed, structural and medical 
imaginJSHSI^ monitoring of pharmaceutical ingredientJ^^, 
and the characterization of material^^i^^. In particu¬ 
lar time-domain terahertz spectroscopy has been proven 
very useful in studying different materia l syst ems at low 
temperatures including superconductor^^^R^, quantum 
magnet^i^®^, excito n sta tes in Ti 02 nanotubed^, and 
topological insulator 

With the photoconductive switches method of TDTS, 
one typically splits an infrared femtosecond laser pulse 
along two paths and sequentially excites a pair of pho¬ 
toconductive “Auston”-switch antennae. An approxi¬ 
mately 1 ps long THz range pulse is emitted by one an¬ 
tenna, focused by mirrors or lenses, transmitted through 
the sample under test, refocused, and then measured at 
the other antenna. By varying the path length differ¬ 
ence of the two paths, the electric field of the transmit¬ 
ted pulse is measured in the time domain. Ratioing the 
Fourier transform of the transmission through the sam¬ 
ple to that of a reference gives the frequency dependent 
complex transmission function in a range that typically of 
order 100 GHz - 3 THz. The complex optical parameters 
of interest e.g. complex dielectric function or the com¬ 
plex index of refraction can then be obtained from the 
complex transmission function. However, one of the con- 
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tinning limiting aspects in using this technique is the typ¬ 
ically strong constraints placed on sample configurations. 
Measurements are generally performed in transmission, 
but in order to get adequate signal to noise (particularly 
at low frequencies) samples must be typically large in 
the transverse direction. In most cases sample dimen¬ 
sions must be of order 4 mm across, which is well in 
excess of the estimates of the diffraction limit of order 
0.6 mm. A principle source of the larger beam spots is 
from the parabolic aberration that originates from finite 
size sources and the off-axis parabolic mirrors which are 
used in most TDTS setup^. 

Because many interesting materials are only found in 
small crystals, it is important to investigate new ways 
of measuring small samples. In this paper we demon¬ 
strate a method to compensate for the large beam sizes 
that come from parabolic aberration and other sources in 
typical terahertz spectrometers. We use a nested coaxial 
parabolic reflector arrangement to couple light through 
a small aperture. The nested coaxial parabolic reflector 
arrangement is designed to have a one-to-one mapping 
of light that takes incoming parallel rays, couples them 
through an aperture as small as 2 mm, and back out 
into parallel rays (Fig. [^. It greatly increases through¬ 
put and signal to noise by a factor of almost 15 over a 
simple aperture, making reliable measurements of sam¬ 
ples as small as 2 mm possible. The reflector arrange¬ 
ment is compact, compatible with cryogenic arrange¬ 
ments, and can be easily manufactured using standard 
CNC machining methods. It may find application in a 
number of circumstances in TDTS measurements in addi¬ 
tion to measuring small samples, including coupling THz 
radiation into subwavelength waveguides^Sl, THz pulse 
shapin^^, and achieving greater overlap between beam 
sizes in optical-pump THz probe spectroscopie^^. An¬ 
other application of this device could be to significantly 
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FIG. 1. (Color online) (a) Three dimensional cross section 
of the nested concentric parabolic reflector; the inset shows a 
photo of the 2 mm reflector (left) and 3 mm reflector (right) 
(b) Schematic view of the two dimensional cross section of 
the nested concentric parabolic reflector. The reflective sur¬ 
face is formed by rotating the intersecting parabola gi and 
g 2 around their axis. F is the focal point of both parabolas, 
whereas / denotes their focal length. The other dimensions 
of the parabolic reflector are labeled by the length I, the inner 
diameter di, the outer diameter d 2 and the aperture size «. 
The yellow cone represents a collimated light beam traveling 
through the device within a ray optics perspective. 


reduce the focal spot size in compact TDTS systems con¬ 
sisting of only two off-axis parabolic mirrors. 


Fresnel equations to be used to determine optical con¬ 
stants of materials from transmission properties. 

In the design the parabolas gi and g 2 are given by the 
formula, 


y = - f)- ( 1 ) 

In these expressions, / corresponds to the focal length 
of the parabolas. The dimensions of the reflector can 
be labeled by several parameters as illustrated in Fig. 
(b). The values chosen for the diameter at the entrance 
of the reflector di and the effective aperture size in the 
middle of the reflector z fix the values for the focal length 
/ and the length of the reflector 1. These parameters are 
respectively given by the following relations. 


2 = 4 / 
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A particular aspect of this design is that - within a ray 
optics perspective - the relative path length difference of 
the reflected rays as compared to the rays going straight 
through the center aperture is a constant independent of 
the ray trajectory. Using Eq. [T} the distance L that a 
light ray travels after hitting the reflective surface until 
it reaches the focal point is calculated as 


II. METHODS 

Fig. g shows the design of the nested concentric 
parabolic reflector. The reflector consists of a cylinder 
containing two nested paraboloids of revolution, facing 
in opposite directions, which are aligned in a coaxial 
manner. The two intersecting paraboloids form a reflect¬ 
ing surface. Fig. (a) and (b) show respectively three 
and a two dimensional cross sections of the reflector. 
The paraboloids of revolution are obtained by rotating 
parabolas gi and g 2 around their axis. The paraboloids 
gi and 52 share the focal point A, as can be seen in Fig. 

(b). Within a ray optics perspective, a collimated light 
beam traveling parallel to the cylindrical surface that en¬ 
ters the top of the reflector, gets reflected by the reflective 
surface of parabola gi and consequently becomes focused 
in the focal point F. Subsequently, the light rays hit the 
reflective surface of parabola 52 and exit the parabolic 
reflector as a one-to-one mapping of the entering beam. 
As this mapping can be motivated on the basis of incom¬ 
ing and outgoing waves being time reversed paths of each 
other, one expects that it is valid within a Gaussian op¬ 
tics perspective as well. Despite the fact that the beam 
is focused, we expect (in a fashion similar to a Gaussian 
beam) that the phase fronts of the propagating beam will 
be flat at the waist e.g in the aperture, allowing general 


L = \/a :2 -P y 2 = V 4/j/ + 4/2 -h =y + 2f (4) 

Here x is the horizontal distance and y is the vertical dis¬ 
tance between the point where the light hits the parabolic 
reflector and the focal point. It follows that the differ¬ 
ence in optical path between two initially parallel rays 
where one is transmitted directly through the reflector 
and the other undergoes two reflections is independent 
of the path through the reflector and given by the ex¬ 
pression 

2(L - y) = 4/ (5) 

Hence, the total path length difference between the re¬ 
flected rays and the light that goes straight through the 
device is determined by the focal length of parabolas and 
interestingly is equivalent exactly to the aperture size z 
itself. 

Two nested concentric parabolic reflectors with differ¬ 
ent dimensions were machined and tested. The outer 
diameter di and inner diameter c ?2 were equal to 10 mm 
and 11 mm respectively. The “2 mm reflector” had an 
aperture size z of 2.0 mm, which dictates a focal length 
of / = 0.5 mm and a reflector length of I = 2.40 cm. 
The “3 mm reflector” had an aperture z of 3.0 mm and 
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FIG. 2. (Color online) (a) Comparison of the TDTS time 
traces of a 12 mm aperture, 2 mm aperture and 2 mm re¬ 
flector, (b) Comparison of the TDTS time traces of a 12 mm 
aperture, 3 mm aperture and 3 mm reflector. 


a focal length / and a reflector length I of 0.75 mm and 
1.52 cm respectively. The nested concentric parabolic 
reflectors were machined from 6061 aluminium alloy us¬ 
ing standard CNC lathe machining methods to an 0.4 
/im surface roughness. The inset in Fig. (a) shows a 
photo of the 2 mm reflector (left) and the 3 mm reflec¬ 
tor (right). In principle our design could be suitable for 
even smaller aperture configurations (<1 mm), but these 
could not be produced on the CNC lathe due to machin¬ 
ing constraints. Sinker-type electrical discharge machin¬ 
ing (EDM) with electropolishing may be useful in realiz¬ 
ing such geometries. The nested concentric parabolic re¬ 
flectors were tested using a standard TDTS system with 
Auston switches and off-axis parabolic mirrors (CAP) in 
an 8f geometry. The parabolic reflector was placed be¬ 
tween OAPs 2 and 3. 


III. RESULTS AND DISCUSSION 

The electric field as a function of time transmitted 
through a 12 mm aperture, a 2 mm aperture and the 
2 mm reflector are shown in Fig. (a). The 12 mm 
aperture time trace contains a single peak, as is expected 


from a typical TDTS time trace. The shape of the 2 mm 
aperture signal is very similar, expect that the maximum 
value of the electric field is 7.3 times smaller than for the 
12 mm aperture. Thus, the focal spot size of the tera¬ 
hertz beam is significantly larger than 2 mm. The time 
trace is more complicated for transmission through the 
reflector. The 2 mm reflector transmission shows a first 
peak at the same location in time as the 2 mm aperture 
that is extremely similar to the 2 mm aperture’s shape, 
except that the maximum electric field of the 2 mm reflec¬ 
tor peak is slightly higher by a factor of 1.3. This effect 
is probably caused by small differences in the size of the 
aperture at the center of the reflector. The initial feature 
is presumably due to the straight transmission through 
the aperture with no reflection off the mirror surfaces. 
Most striking in the reflector spectrum however is the 
second feature that shows up slightly later in time. The 
shape of this second feature is somewhat similar to the 
shape of the first peak, but the second peak is broader 
in time, larger in magnitude and phase shifted as com¬ 
pared to the first peak. The time difference between the 
first and the second peaks is 7.0 ps, which corresponds 
to a path length increase of 2.1 mm. This path length 
difference is approximately equal to the aperture size z 
as predicted in Eq. above for the reflected peak. It can 
be concluded that the second peak is due to radiation be¬ 
ing reflected by the parabolic surface and subsequently 
being passed through the aperture at the center of the re¬ 
flector. The maximum electric field of the second peak, 
which will be called the reflector peak hereafter, is 2.6 
times as large as the first peak. 

Fig. i (b) shows the electric field of the transmission 
through respectively a 12 mm aperture, a 3 mm aperture 
and the 3 mm reflector as a function of time. Similar 
trends as for the 2 mm reflector can be observed and 
in this case the time difference between the first peak 
and the reflector peak is larger and corresponds to 3.08 
mm. Again, this is approximately equal to the size of the 
aperture z as has been predicted above. The maximum 
electric field of the second peak is now 1.4 times larger 
than the first peak. The smaller ratio between these two 
peaks in the case of the 3 mm reflector indicates that 
this reflector is less effective at coupling additional light 
through the aperture in the center. This is presumably 
expected because a larger percentage of the signal can 
make it through on the first pass with no reflection. We 
note that additional polishing of the 3 mm reflector to 
a mirror finish did not increase the throughput, but in¬ 
stead induced an additional and erratic phase shift in 
the signal. This is probably caused by deviations in the 
parabolic shape which were induced in the polishing pro¬ 
cess. 

Fig.0 shows a comparison of the spectral amplitudes, 
found by taking the Fourier Transform of the time traces 
in Fig. of the 12 mm aperture, 2 mm aperture, and 2 
and 3 mm reflectors. It can be seen that the 12 mm aper¬ 
ture signal has a much higher spectral amplitude than 
the 2 and 3 mm aperture. This difference in magnitude 
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FIG. 3. (Color online) (a) Comparison of the spectral ampli¬ 
tudes of the throughput of a 12 mm aperture, 2 mm reflector 
and 2 mm parabolic reflector, (b) comparison of the spectral 
amplitudes of the throughput of a 12 mm aperture, 3 mm 
reflector and 3 mm parabolic reflector. The reflector spec¬ 
tral amplitude is analyzed by taking Fourier transform of the 
whole signal (“both peaks”) or only of the second peak (“only 
reflector peak”). 

is most pronounced at low frequencies (^0.1 up to ~0.8 
THz), because of the large focal spot size of the low fre¬ 
quency light. The focal spot size of the high frequency 
light is smaller, therefore the difference in spectral mag¬ 
nitude between the apertures is smaller. The reflector 
signals were analyzed in two different ways. The spec¬ 
tral amplitude obtained by taking the Fourier transform 
of the complete reflector signal shows strong oscillations. 
These oscillations are caused by interference between the 
first peak and the reflector peak in the time trace. These 
oscillations can be removed by only analyzing the reflec¬ 
tor peak. This has been achieved by subtracting the 
rescaled aperture time trace from the reflector time trace 
before taking the Fourier transform. As can be seen in 
Fig. [3 this removes most of the oscillations. The os¬ 
cillations cannot be completely removed in this fashion 
because the shape of the first peak in the reflector and 
aperture time trace are not identical. 

The maximum spectral amplitude in the case of the 
reflector is 15 times the amplitude of a regular aperture 
of the same size in the case of the 2 mm reflector and 5 


times for the 3 mm reflector at 0.27 THz. Between 0.20 
and 0.35 THz the spectral amplitude of the reflector is of 
the same order as the the 12 mm aperture. The reflectors 
significantly increase the signal up to a frequency of 0.7 
THz. The reflector peak is broader in time than the first 
peak as can be seen in Fig. [3 A peak which is broader in 
time is composed of lower frequencies and this is indeed 
the case as can be seen in the spectrum. Strictly speaking 
this low frequency radiation does not obey the traditional 
laws of ray optics, but instead forms a Gaussian beam 
waist with a large spot size when entering the reflector. 
The reflectors couple this low frequency light through 
the hole very efficiently. Light with a frequency above 
~1 THz has a much smaller spot size. Most light in this 
domain makes it through the aperture without reflecting. 
Therefore, the reflector does not significantly increase the 
spectral amplitude in this region. 


IV. CONCLUSION 

We have demonstrated an effective way of increasing 
the light throughput of a small aperture. It can be con¬ 
cluded that the novel nested concentric parabolic reflec¬ 
tor significantly reduces the effective focal spot size of 
low frequency terahertz radiation by coupling it through 
a small aperture. The results suggest that a reflector 
with a smaller aperture than demonstrated in the paper 
might be even more efficient at increasing the throughput 
of light. 

This technique could enable the study of smaller sam¬ 
ples than has been possible so far by means of TDTS 
experiments, by placing a sample in a slot machined in 
the middle of the parabolic reflector. In order to main¬ 
tain the one-to-one mapping of in-going to out-going rays 
and the proper focussing at least approximately, the de¬ 
sign would presumably have to account for the additional 
optical path length in the sample under test. Moreover in 
order to be able to analyze the data, it is important that 
the reflector peak does not overlap with the reflections 
from Fabry-Perot interference from multiple reflections 
in a substrate. Thus, optically thick samples and/or sub¬ 
strates should be used to clearly separate these peaks. 
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